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The spin dynamics of the semiclassical Heisenberg model with uniaxial anisotropy, on 
the layered triangular lattice with antiferromagnetic coupling for both intralayer nearest 
neighbor interaction and interlayer interaction is studied both in the ordered phase and in the 
paramagnetic phase, using the Monte Carlo-molecular dynamics technique. The important 
quantities calculated are the full dynamic structure function S{c[,uj), the chiral dynamic 
structure function Schii^^), the static order parameter and some thermodynamic quantities. 
Our results show the existence of propagating modes corresponding to both ^(q, a;) and 
Schii^^) in the ordered phase, supporting the recent conjectures. Our results for the static 
properties show the magnetic ordering in each layer to be of coplanar 3-sublattice type 
deviating from 120 degree structure. In the presence of magnetic trimerization however, 
we find the 3-sublattice structure to be weakened alongwith the tendency towards non- 
coplanarity of the spins. 

1 Introduction 

The layered triangular antiferromagnets with crystal field induced single ion uniaxial anisotropy 
are very important both from the viewpoint of behavior of magnetic materials as well as from 
the interest in spin models and universality classes in statistical mechanics^'^. The materials 
CsMnBr2, VCI2, VBr2 are some of the materials belonging to this class. They exhibit a 
complex competition involving geometric frustration and spin component degeneracy along 
the easy axis. The neutron scattering experiments on these materials both in the ordered 
phase as well as in the disordered phase by the experimentalists^ provide a new impetus 
for theoretical study. The physical quantities of principal interest are the nature of static 
magnetic ordering and spin dynamics. In particular there are very interesting questions 
regarding the coplanarity and the sublattice structure of the magnetic ordering^'^. In dy- 
namics there are speculations about the chiral mode'^. The only reliable way at present 
to calculate the static correlations of the Heisenberg model is by the technique of Monte 



Carlo (MC) simulation. Furthermore the spin dynamics of this model can also be accurately 
studied only by using the technique of MC combined with Molecular Dynamics (MD). We 
have been motivated by the recent experiments to compute dynamic spin chirality correlator 
structure factor Schii^^) and total dynamic structure factor S{q, uj) for this model appropriate 
to layered triangular antiferromagnets both below and above Neel temperature. Besides, we 
also analyse the detailed nature of the magnetization vector(the conventional order param- 
eter) and the effect of magnetic trimerization on the ordering pattern. Here we carry out 
our MCMD calculation with model parameters appropriate to Mn based systems with spin 
value(S) equal to 2. We apply quantum correction to our semiclassical MCMD calculation 
to make our theoretical results comparable with the experimental ones. 

The plan of the paper is as follows:- In Sec. II we describe the model and the quantities 
of interest which we calculate; in Sec. Ill we very briefly describe the MCMD technique 
and our calculational procedure; in Section IV we display all the important results from our 
calculation; and finally in Sec. V we discuss some implications of our results. 

2 Model and Mathematical formulation 

The model appropriate for our investigation is 

y: j^jS,.sj+Dj:sf (1) 

where Jij is positive and takes the value Ji for ij being nearest neighbors (i.e. for 
intralayer coupling) and the value J2 for ij being next nearest neighbors (i.e. interlayer 
couphng); D is -ve. The definitions of most of the quantities of interest viz. S{q,uj), M etc. 
are very common and can be found in Ref.5. The new and very interesting quantity studied 
here viz. Schii^^) is defined as :- 

00 

Sehiioo) = J dtC,hi{t)e-"^' (2) 
—00 

where Cchi{t) is the scalar spin chirality correlator. It is defined as :- 

CeMit) = (x(O)x(t)), X(t) = E X^d^), X(0) = E X'tmniO) (3) 

(ijk) {Imn) 

where {ijk) and {Imn) represent the plaquettes in the triangular lattice. The total lattice 
scalar spin chirality x(^) has been defined in terms of the local scalar spin chirality X^jki^)- 
The ( ) occurring in the equation for Cchi{t) have been used to represent the thermodynamic 
average. The sum in the definition of lattice chirality is taken over a plaquette in the anti- 
clockwise direction. The local spin chirality itself is defined as:- 

X^.fc = {SrS.XS,) (4) 

The other static quantity introduced in this paper is inverse degree of coplanarity 7 defined 
as:- 

7 = [Ma-MbXMc]/[\Ma\\Mb\\Mc\] (5) 

where M^,]VIs and M^; are the magnetization vector order parameters in the three sublat- 
tices A,B and C respectively. We evaluate 7 for our model both in the absence and in the 
presence of trimerization. 
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3 The MCMD Approach and Calculations 



The basic equation governing the spin dynamics is the following: 



where the square bracket denotes the commutator. This equation leads to more simplified 
equations of motion :- 

^ - f E^^.(s.^s,) + [-x^^. + y^^.] (7) 

The more details regarding these can be found in Ref.s 5 and 6. Wc take the semiclassical 
version of the above equation for spin dynamics^. We then apply quantum correction to 
evaluated semiclassical Scui'^) and S{q,ui) following Windsor's prescription^. The corrected 
Schii'^) takes the following form in this prescription:- 

STiu^) = [-—^]ScHf\u;) (8) 

1 + e 2,r 

We use the Metropolis mcthod^'^ for MC. We adopt the sequential updating procedure 
for MC updating in our calculation. As our model is highly frustrated, it takes long time 
to equilibriate. Therefore we discard all the spin configurations with MC age below 5000 
MC steps/spin (MCS/spin) for static calculations involving thermodynamic averages. Wc 
store spin configurations with MC ages between 8000 — 10, 000 MCS/spin, 14, 000 — 16, 000 
MCS/spin and 40,000 — 42,000 MCS/spin for various calculations and investigations, de- 
pending upon the temperatures of interest and the model parameters. In particular these 
configurations are also used as boundary conditions for MD calculations. The configurations 
used for static calculations numbered 1000 and those for MD numbered 5. Wc use the lattice 
of size 9X9X9 containing 729 spins. The magnetic ordering temperature (T/v) was deter- 
mined by locating the peak for the specific heat Cy{T) calculated by MC. We had first carried 
out our MC calculations with parameters corresponding to those of Kawamura^ and found 
a very good agreement between our results and his. This was done as a reliability check on 
our calculations. The equation of motion in our model in the semiclassical approximation 
takes the form :- 

at ., J\ I \ 

jynn) j(nnn) 

where to is the natural time unit in our problem, where 

to-'^^MSis + i))-^ (10) 

The vector Si represents a unit spin vector at the lattice site i with the operator replaced 
by the classical vector Si{S{S + l))^ ■ We calculate Schiiy^) as a function of uj for our model 
at different temperatures both in the ordered phase as well as in the disordered phase. Later 
we also carried out MCMD calculation in the presence of a small magnetic field of strength 
5 Tesla applied along x direction to verify the intrinsicity of this chiral mode. 
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Regarding the study on trimerization, wc take the foUowing model, 

= Ji(l + (-l)*A) 



(11) 



where 5 takes the value even for the plaquette involving the points [{x,y),{x + l,y + 
1), {x,y + 1)] and the value odd otherwise; A is a positive quantity denoting the enhance- 
ment/suppression factor of the nearest neighbor (nn) coupling and jp"™-*^ is the modified nn 
coupling in the presence of trimerization. Moreover there are some special sets of values for 
X and y where this trimerization occurs^. We adopt the same pattern in our calculation. All 
the simulation work reported in this paper was carried out using two machines. The MC 
calculations were done on Dec Alpha (IISC) and the MD calculations were carried out on a 
networked Linux machine at Bose Centre. The total duration of each dynamical run in MD 
denoted by tmax was chosen as 75to to get resolution width of 0.8 mev. 

4 Results 

In this section we present all our important results, for both the static and the dynamic 
properties. We have done calculations with two sets of model parameters which have same 
values for Ji and J2 but differ in the values of D. We calculated for both this model by 
MC calculation. The two sets of model parameters used by us are:- 

Ji = 57.0K, J2 = 5.7K, D = -5.7K (12) 

for the first system and 

Ji = 57.0K, J2 = 5.7K, D = -8.55X (13) 

for the second system. From our MC calculations we determine Tjv for both the systems. 
They are found to be:- 

= 66K, = 68K (14) 
For trimerization study we chose the system parameters as 

JN = 1.0, JNN = 1.0, D = 0.0 (15) 

This gives 

Tn = 1.2 (16) 

for the untrimerizcd case. For the trimerized case we choose A to be equal to 0.2. This gives 
the value of T/v to be close to 1.0. 

The sublattice magnetization vectors for all the 9 layers in the ordered phase are deter- 
mined from our MC calculations at the temperature of 20K. These calculations have been 
done for 2 MC samples prepared independently with MC ages varying between 14, 000 and 
16,000 MCS/spin. Wc have done this study only for the first system. The results for total 
magnetization vectors (taking into account the all 3 sublattice magnetization vectors) for 
sample 1 and sample 2 are presented in tables 1 and 2 respectively. It should be pointed 
out that the true values of the magnitudes of the magnetization can be obtained by dividing 
the corresponding magnitudes occurring in the tables 1 and 2 by a factor of 3. We have 
also studied the thermal evolution of 7 and it is presented in table 3. Moreover the effect of 
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Table 1: Total magnetization vectors for 9 layers for MC sample 1 corresponding to system 
1 at 20K; M denotes the total magnetization, M^y the total magnetization vector in the xy 
plane etc., (f)xy denotes the angle of the vector M^y with the x-axis, Vm represents the 3-d 
volume (with the sign) made by the 3 components of M and the angle has been given 
in units of tt. 
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Figure 1: Schi{i-^) vs. u corresponding to system 1 at temperature of 0.3T^; for layer 1 (left 
graph) and for averaging over all layers (right graph) 



trimerization on 7 has also been studied for our model with D put to zero. This result is 
displayed in table 4. 

Regarding the dynamical properties, we present both Schii^^) vs. u and S{q,uj) vs. u 
for constant q-scans, obtained from our MCMD calculations, in the form of various figures. 
We consider both individual layers as well as layer averaged value in these graphs. We have 
carried out these calculations for our model with parameters appropriate to both first system 
and second system. The q vectors have been taken along both < 110 > as well as < 111 > 
direction for our calculation of S'(q, a;). 



5 Discussion 

From tables 1 and 2 we see that there is a non-trivial correlations amongst M occurring in 
the 9 layers. Firstly shows anti-parallel correlations between adjacent layers, excepting 
two layers, for both the MC samples. The same thing cannot however be said regarding M^y. 
In sample 2 it mostly shows an antiparallel arrangement but not so for sample 1. Moreover 
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Table 2: Same as table 1 but for MC sample 2 
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Table 3: 7(7") for system 1 in the untrimerized case 
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Table 4: 7(T) for system 1 in the trimerized case 



Temperature (T) 
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Table 5: < ujp >iayav (T) vs. T for both system 1 and system 2 
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Figure 2: Schi{^) vs. uj corresponding to system 2; for layer 1 at temperature of 0.75TAr 
(left most graph), for layer 3 at temperature of 0.75Tn (centre graph) and for any layer at 
temperature of 0.25r;v (right most graph) 




Figure 3: Schii^) vs. ou corresponding to system 1; for layer 1 in the presence of magnetic 
field of 5 Tesla at temperature of 0.5T/v (left graph) and for any layer in the absence of any 
magnetic field at temperature of 0.3T/v (right graph) 
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Figure 4: S{q,uj) vs. uj corresponding to system 1 for layer 1 at temperature of 0.3T/v; 
values of q are (0, 0, 0) (left graph) and (6, 6, 0) (right graph) 



the 



in both sample 1 as well as in sample 2, a few layers exhibit a large suppression of total 
magnetic moment. The value of 7 also shows variation with layer index hut indicates a "fairly 
high degree of coplanarity" for the strictly non-coplanar 3-sublattice structure of magnetic 
ordering, for all the layers for both the MC samples. This result is of much significance in the 
presence of uniaxial anisotropy. In contrast to the case with pure Heisenberg model however, 
the angles between the 3 sublattice ordering directions are not 120 degrees in general in our 
case. Prom the results presented in tabic 3 wc sec that the degree of coplanarity (layer 
averaged) for system 1, has an oscillatory behavior as the temperature is lowered from T^. 
From table 4 we see that trimerization causes coplanarity to decrease a bit in the absence 
of the anisotropic term. Moreover we also observed that 3-sublattice structure is much 
weaker i.e. much less distinct in the case of trimerized model compared to the untrimerized 
case. These results support the conjectures of the experimentalists^. Wc also sec that with 
decrease in temperature the degree of coplanarity increases even for the trimerized system. 
This study also brought out the fact that trimerization causes T/v to drop. 

Regarding spin dynamics, we see many more interesting features in Schi{<^) than in 
5'(q, a;). There exists a well defined propagating mode for Schi{(^) in the ordered phase 
for both the sets of system parameters, as are found from the figures 1 and 2. These graphs 
describe the behavior of Schi{'-o) as a function of uo for some particular layers. The frequencies 
of propgating modes i.e. the values of the frequencies at which the curves show the peaks, 
are found to be dependent on the layers, for both the systems. The layer averaged peak fre- 
quencies have also been calculated. These frequencies denoted by < 0;^ > have been studied 
as a ftmction of temperature for both the systems. The results of this study arc displayed in 
table 5. These modes disappear in the disordered phase i.e. for T above Tjv, excepting in the 
case of system 2. It should be pointed out that these < uOp >'s arc in general different from 
the peak position of the layer averaged Schi{^)- Another notable feature in these figures is 
the existence of a dip at higher value of uj and that occurs in the regime of -ve values for 
Schi- This regime probably describes the spin fiipping process causing the chirality to be 
reversed at high energies. From the analysis of these results we see a very clear thermal 
evolution of the chiral mode for system 1. The softening of averaged peak frequency with 
increase in temperature is seen quite distinctly for systeml. For system 2 with higher value 
for the magnitude of D however, we hardly see any thermal evolution for the peak frequency. 
This is probably due to the freezing of spin dynamics by the rather large value of D in the 
low u region. It is to be noted that our Schi{i^) is actually SchiiQ = 0,cj). The intrinsicity 
i.e. non-spurious nature of these modes were verified by applying a magnetic field, as has 
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been described earlier. The peaks were found to essentially vanish and become "shoulders" 
in the presence of applied field, in contrast to the zero field situation. This is brought out 
in figure 3. The existence of these chiral modes support the theoretical conjecture'^. At low 
temperature the layer averaged peak frequency corresponding to this chiral mode, seems to 
scale as D*^'^, with the same values for J's and 7^. It is important to add that these chiral 
modes disappear in the disordered phase i.e. for T above T^, as we observed in our MCMD 
calculation for system 1 at the temperature of l.lTjv. We however have not presented the 
results here. 

The total spin dynamics described by the full dynamic structure function ^(q, u) as a 
function of uj in the constant-q scan, exhibit spin wave like collective modes in the ordered 
phase for both the systems. These modes however look quite dispersionless, both with 
respect to the magnitude and direction of q. They do not show any appreciable temperature 
dependence (softening) in the ordered phase either. Besides this finite frequency peak, we 
also see a central peak i.e. a peak at zero frequency. The relative intensities at these two 
peaks are however a function of temperature and the value of q. These results are presented 
in the figure 4. The finite frequency peak is quite sharp at low temperature and is of much 
higher intensity than the central peak for large values of |gf|, as we can see from these figures. 
The central peak arises probably from the very slow dynamics of the longitudinal components 
of the spin, the longitudinal directions being the 3-sublattice ordering directions. For low 
values of |g| however, the central peak dominates over the finite 00 peak in terms of intensity. 
The apparent dispersionless and gapped behavior of these "spin waves" are probably due to 
the presence of D term. The observed insensitivity of these modes to temperature in the 
ordered phase also arises from this anisotropy term, most likely. 

It is quite striking to observe the contrast in the behavior of Schiii-^) and S'(q, a;) partic- 
ularly with respect to thermal dependence. This probably refiects the lesser influence the 
single ion anisotropy has on 2-plaquette correlations than on 2-spin correlations. Moreover 
no central peak occurs in the presence of the finite frequency peak for Schi{i^) in contrast 
to the case with 5'(q, a;). It is to be noted that the position of the "spin wave" peak looks 
consistent with the magnitude of Ji. It is extremely interesting to notice that the position 
of the finite frequency peak for S'(q, cj) is very close to the position of the dip for Schi{^^)- 
On the other hand the peak position of Schiiy^) is close to the positon of the central peak for 
S'(q, cu). This probably implies that the dominant contribution to the total spin dynamics 
here comes from the chiral dynamics. It would also be interesting to study the dynamic 
structure factor corresponding to the staggered chiral correlator Schi{f\,^) by MCMD. In 
brief the spin dynamics of a geometrically frustrated system in the presence of single ion 
anisotropy is very rich. It would be quite exciting to determine the role and importance 
of quantum fluctuations in these systems {S — 2) by comparing our semiclassical results 
for dynamic structure factors with those from inelastic neutron scattering experiments in a 
detailed quantitative way . Unfortunately detailed experimental results for these systems 
are still awaited. 
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